The  Astronomical  Journal,  142:190  (13pp),  2011  December 

©  2011.  The  American  Astronomical  Society.  All  rights  reserved.  Printed  in  the  U.S.A. 


doi:  10.1088/0004-6256/142/6/190 


EXPLORING  THE  VARIABLE  SKY  WITH  LINEAR.  I.  PHOTOMETRIC 
RECALIBRATION  WITH  THE  SLOAN  DIGITAL  SKY  SURVEY 

Branimir  Sesar1,2,  J.  Scott  Stuart3,  Zeljko  Ivezic2,  Dylan  P.  Morgan2,4, 

Andrew  C.  Becker2,  and  Przemyslaw  Wozniak5 

1  Division  of  Physics,  Mathematics  and  Astronomy,  Caltech,  Pasadena,  CA  91125,  USA 

2  Department  of  Astronomy,  University  of  Washington,  Seattle,  WA  98195-1580,  USA 

3  Lincoln  Laboratory,  Massachusetts  Institute  of  Technology,  Lexington,  MA  02420-9108,  USA 
4  Department  of  Astronomy,  Boston  University,  Boston,  MA  02139,  USA 
5  Los  Alamos  National  Laboratory,  Los  Alamos,  NM  87545-0001,  USA 
Received  2011  July  13;  accepted  2011  September  23;  published  2011  November  4 

ABSTRACT 

We  describe  photometric  recalibration  of  data  obtained  by  the  asteroid  survey  LINEAR.  Although  LINEAR  was 
designed  for  astrometric  discovery  of  moving  objects,  the  data  set  described  here  contains  over  5  billion  photometric 
measurements  for  about  25  million  objects,  mostly  stars.  We  use  Sloan  Digital  Sky  Survey  (SDSS)  data  from  the 
overlapping  ~  10,000  deg2  of  sky  to  recalibrate  LINEAR  photometry  and  achieve  errors  of  0.03  mag  for  sources 
not  limited  by  photon  statistics  with  errors  of  0.2  mag  at  r  ~  18.  With  its  200  observations  per  object  on  average, 
LINEAR  data  provide  time  domain  information  for  the  brightest  four  magnitudes  of  the  SDSS  survey.  At  the  same 
time,  LINEAR  extends  the  deepest  similar  wide-area  variability  survey,  the  Northern  Sky  Variability  Survey,  by 
3  mag.  We  briefly  discuss  the  properties  of  about  7000  visually  confirmed  periodic  variables,  dominated  by  roughly 
equal  fractions  of  RR  Lyrae  stars  and  eclipsing  binary  stars,  and  analyze  their  distribution  in  optical  and  infrared 
color-color  diagrams.  The  LINEAR  data  set  is  publicly  available  from  the  SkyDOT  Web  site. 
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1.  INTRODUCTION 

Variability  is  an  important  phenomenon  in  astrophysical 
studies  of  structure  and  evolution,  both  stellar  and  Galactic 
and  extragalactic.  Some  variable  stars,  such  as  RR  Lyrae 
stars,  are  an  excellent  tool  for  studying  the  Galaxy.  Being 
nearly  standard  candles  (thus  making  distance  determination 
relatively  straightforward)  and  being  intrinsically  bright,  they 
are  a  particularly  suitable  tracer  of  Galactic  structure  (e.g.,  Sesar 
et  al.  2007,  and  references  therein).  Similarly,  eclipsing  binaries 
can  be  used  as  distance  indicators  (Guinan  et  al.  1998),  are 
excellent  probes  of  stellar  physics  (Prsa  &  Zwitter  2005),  and 
offer  a  unique  method  for  measuring  stellar  masses  and  other 
parameters  as  well  (e.g.,  Andersen  1991;  Torres  et  al.  2010). 

Despite  the  importance  of  variability,  the  variable  optical  sky 
remains  largely  unexplored  and  poorly  quantified,  especially 
at  the  faint  end  ( V  >  15).  To  what  degree  do  different 
variable  populations  contribute  to  the  overall  variability,  how 
are  they  distributed  in  magnitude  and  color,  and  what  are 
the  characteristic  timescales  and  dominant  mechanisms  of 
variability  are  just  some  of  the  questions  that  still  remain 
unanswered.  To  address  these  questions,  several  contemporary 
projects  aimed  at  regular  monitoring  of  the  optical  sky  were 
started.  The  four  most  prominent  wide-area  (more  than  a  few 
thousand  deg2)  variability  surveys  in  terms  of  depth  and  cadence 
are  as  follows. 

1 .  ROTSE-I  (Akerlof  et  al.  2000)  monitored  the  entire  observ¬ 
able  sky  twice  a  night  from  V  =  10  to  a  limit  of  V  =  15.5. 
The  Northern  Sky  Variability  Survey  (NS VS;  Wozniak  et  al. 
2004)  is  based  on  ROTSE-I  data  and  contains  light  curves 
for  about  14  million  objects,  with  100-500  measurements 
per  object  collected  over  one  year.  Typical  photometric  pre¬ 
cision  of  this  data  set  is  about  0.02  mag  for  sources  not 
limited  by  photon  statistics. 


2.  The  All  Sky  Automated  Survey  (AS AS;  Pojmanski  2002) 
monitors  the  entire  southern  and  part  of  the  northern  sky 
(8  <  +28°)  to  a  limit  of  V  =  14,  including  about  15  million 
stars.  Typical  photometric  precision  is  about  0.05  mag.  The 
third  release  of  the  AS  AS  Catalog  of  Variable  Stars  contains 
close  to  50,000  variable  stars.  About  80%  of  these  are  new 
discoveries  (not  previously  cataloged). 

3.  The  Lowell  Observatory  Near  Earth  Objects  Survey  Phase 
I  (LONEOS-I;  Miceli  et  al.  2008)  provides  photometric 
data  for  1430  deg2  of  northern  sky  that  has  been  imaged  at 
least  28  times  between  1998  and  2000.  The  LONEOS-I 
camera  used  no  bandpass  filter  and  reached  a  depth  of 
R  ~  18.5.  Typical  photometric  precision  of  this  data  set  is 
about  0.02  mag  for  sources  not  limited  by  photon  statistics. 

4.  The  Palomar  Transient  Factory  (PTF;  Law  et  al.  2009 
and  Rau  et  al.  2009)  is  an  ongoing  wide-area,  two-band 
(SDSS-g'  and  Mould-7?  filters),  deep  (R  ~  20.6,  g'  ~  21.3) 
survey  aimed  at  systematic  exploration  of  the  optical 
transient  sky.  As  of  fall  2011,  PTF  has  observed  about 
7300  deg2  of  northern  sky  at  least  30  times  in  the  Mould-7? 
band  (~1800  deg2  of  sky  with  more  than  100  epochs). 
Typical  photometric  precision  of  this  data  set  is  better  than 
0.01  for  sources  not  limited  by  photon  statistics. 

A  comprehensive  review  of  past  and  ongoing  variability 
surveys  can  be  found  in  Becker  et  al.  (2004). 

With  hundreds  of  observations  per  object,  the  NS  VS  and 
AS  AS  data  sets  represent  unique  resources.  Although  a  number 
of  massive  synoptic  surveys,  such  as  Pan-STARRS  (Kaiser  et  al. 
2002)  and  LSST  (the  Large  Synoptic  Survey  Telescope;  Ivezic 
et  al.  2008),  will  significantly  improve  our  knowledge  of  the 
variable  faint  sky,  they  will  not  obtain  hundreds  of  observations 
per  object  any  time  soon. 
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A  downside  of  NSVS,  ASAS,  LONEOS,  and  PTF  data 
sets  is  the  lack  of  precise  multi-color  photometry,  which  can 
provide  valuable  information  about  sources  in  addition  to  light 
curve  characteristics  (Covey  et  al.  2007).  The  Sloan  Digital 
Sky  Survey  (SDSS;  York  et  al.  2000)  has  obtained  five-band 
photometry  for  hundreds  of  millions  of  stars,  but  its  temporal 
coverage  is  very  poor.  The  only  region  on  the  sky  with  more 
than  10  observations  per  source  is  the  ~300  deg2  large  equatorial 
region  called  Stripe  82  (Sesar  et  al.  2007,  and  references  therein). 
At  the  same  time,  NSVS  and  ASAS  are  not  deep  enough  to  make 
efficient  use  of  SDSS  photometry.  This  goal  to  combine  a  large- 
area  variability  survey  providing  hundreds  of  observations  with 
exquisite  SDSS  photometry  was  recently  made  reachable  by  the 
LINEAR  survey. 

The  MIT  Lincoln  Laboratory  (MITLL)  has  operated  the 
Lincoln  Near-Earth  Asteroid  Research  (LINEAR)  program 
since  1998  to  search  for  asteroids.  The  focus  of  the  program 
is  to  discover  and  track  near-Earth  asteroids  (NEAs)  larger 
than  1  km  in  diameter,  but  LINEAR  has  also  discovered  and 
tracked  many  main-belt  asteroids  and  smaller  NEAs.  Since  late 
2002,  the  LINEAR  program  has  maintained  an  archive  of  all 
of  its  imagery  data.  This  archive  now  contains  approximately 
6  million  images  of  the  sky,  most  of  which  are  5  megapixel 
images  covering  2  deg2.  The  LINEAR  image  archive  contains 
a  unique  combination  of  depth  ( V  <  18),  coverage  area,  and 
cadence  that  is  potentially  useful  in  exploring  a  number  of  time- 
varying  astronomical  phenomena. 

Since  the  LINEAR  project  was  focused  on  astrometric  survey, 
standard  stars  required  for  precise  photometric  calibration  were 
often  unavailable.  Furthermore,  a  significant  fraction  of  the  data 
was  taken  in  non-photometric  conditions.  The  main  purpose  of 
this  paper  is  to  describe  photometric  recalibration  of  LINEAR 
data  using  stars  observed  by  SDSS  as  photometric  standards. 
Companion  papers  are  discussed  in  Section  6. 

In  Section  2  we  describe  the  LINEAR  survey  in  more  detail 
and  discuss  recalibration  procedures  in  Section  3.  The  recali¬ 
brated  data  set  is  characterized  in  Section  4  and  a  preliminary 
analysis  of  about  7000  confirmed  periodic  variables  is  presented 
in  Section  5.  Our  results  are  summarized  in  Section  6. 

2.  LINEAR  SURVEY 
2.7.  Hardware  Characteristics 

The  LINEAR  program  operates  two  telescopes  at  the  Ex¬ 
perimental  Test  Site  located  within  the  US  Army  White  Sands 
Missile  Range  in  central  New  Mexico  at  an  altitude  of  1506  m 
above  sea  level  and  a  latitude  of  33?  8N.  The  program  uses  two  es¬ 
sentially  identical  equatorial-mounted  telescopes  of  the  Ground- 
based  Electro-Optical  Deep  Space  Surveillance  (GEODSS) 
type,  usually  denoted  LI  and  L2.  Each  telescope  is  a  folded 
prime-focus  design  with  a  primary  mirror  diameter  of  1.016  m 
and  //#  of  2.15  (Stokes  et  al.  2000). 

The  two  telescopes  of  the  LINEAR  project  are  each  equipped 
with  a  CCD  camera  developed  at  MITLL  using  CCDs  de¬ 
veloped  and  manufactured  by  MITLL  for  the  US  Air  Force 
GEODSS  system  (Burke  et  al.  1998).  These  CCID-16  sen¬ 
sors  are  5  megapixel  (2560  x  1960),  back-illuminated,  frame- 
transfer  CCDs  with  fast  readout,  low  dark  current,  and  low  read 
noise.  The  CCID-16  sensors  have  eight  readout  regions  (here¬ 
after  referred  to  as  cells)  and  have  high  quantum  efficiency  over 
a  broad  part  of  the  visible  and  near-infrared  spectrum  (with 
overall  solar- weighted  quantum  efficiency  of  65%  and  peak  ef¬ 
ficiency  of  96%  at  620  nm,  see  Figure  1).  The  cameras  have  no 


Figure  1.  LINEAR  (dashed  line)  and  SDSS  ugriz  spectral  response  curves 
(solid  lines).  The  response  curve  for  LINEAR  includes  the  spectral  quantum 
efficiency  of  the  back-illuminated  CCD,  reflectivity  of  the  aluminum  mirrors, 
and  atmospheric  transmissivity  at  airmass  of  1.3.  The  response  curves  for  SDSS 
include  the  same  atmospheric  absorption. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

spectral  filters.  The  combination  of  the  CCID-16  camera  and 
the  1  m  GEODSS  telescope  produces  an  instantaneous  field  of 
view  of  1?60  x  1?23  (~2  deg2),  with  2"25  pixels.  In  Section  3, 
the  combination  of  CCID-16  sensors  and  LI  and  L2  telescopes 
will  be  referred  to  as  the  LI  and  L2  setups. 

For  a  few  months  in  2003  (approximately  August  to  Novem¬ 
ber)  the  camera  on  the  L2  telescope  was  replaced  with  a  very 
similar  camera  containing  a  smaller  format  sensor  called  the 
CCID-10.  This  smaller  sensor  is  identical  to  the  CCID-16’s  in 
terms  of  manufacture  and  processing,  readout  electronics  and 
settings,  spectral  quantum  efficiency,  and  pixel  size.  The  only 
difference  between  the  smaller  format  camera  and  the  larger  for¬ 
mat  camera  is  the  number  of  pixels,  1024  x  1024  for  the  small 
camera,  and  the  number  of  readout  regions  (4  cells).  During  the 
four  months  in  2003  that  L2  was  operating  with  the  CCID-10 
camera,  the  field  of  view  was  about  0?64  x  0?64  (~0.4  deg2).  In 
Section  3,  the  combination  of  CCID-10  sensor  and  L2  telescope 
will  be  referred  to  as  the  small  setup. 

2.2.  Observing  Conditions  and  Cadence 

The  telescope  site  has  reasonably  dark  skies,  approaching 
liy  =  21  mag  arcsec-2  on  the  best  nights.  However,  the 
LINEAR  program  operates  on  nights  with  suboptimal  condi¬ 
tions  including  haze,  bright  moonlight  (to  within  3  days  of  the 
full  moon),  scattered  clouds,  and  airglow.  LINEAR  images  are 
frequently  taken  under  skies  as  bright  as  {iy  —  16  mag  arcsec-2. 
The  point-spread  function  includes  unknown  contributions  from 
atmospheric  seeing  and  various  system  contributions  (e.g.,  de¬ 
focus  and  dome  airflow)  and  is  typically  around  5"  FWHM, 
resulting  in  slightly  undersampled  images. 

The  image  integration  time  is  allowed  to  vary  over  the  seasons 
so  that  a  roughly  constant  number  of  search  fields  fills  the 
available  dark  time.  During  the  earlier  days  of  the  LINEAR 
program  (roughly  1998  to  2006),  the  integration  time  was 
10-12  s  on  long  winter  nights  and  3-5  s  on  short  summer  nights. 
In  recent  years  (roughly  2007  to  2009),  the  integration  time  was 
increased  to  15-18  s  in  the  winter  and  8-13  s  in  the  summer. 
Under  these  observing  conditions,  the  magnitude  for  which  the 
magnitude  error  is  ~0.2  mag  (the  5a  detection  limit)  is  about 
1 8  mag  on  average. 

As  shown  in  Figure  2,  the  cadence  of  LINEAR  observations 
(time  between  revisits  of  the  same  patch  of  sky)  is  fairly  uniform, 
with  a  peak  at  ~15  minutes  and  a  gap  at  ~8  hr.  A  more 
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log(seconds  between  revisits) 


Figure  2.  Cadence  of  LINEAR  observations  measured  as  the  time  between 
revisits  of  the  same  patch  of  sky.  From  left  to  right,  the  dashed  lines  show 
1  minute,  1  hr,  1  day,  1  week,  1  month,  and  1  year  cadences.  The  LINEAR 
temporal  coverage  is  fairly  uniform,  with  a  peak  at  ~1000  s  corresponding  to 
the  main  15  minute  cadence. 

detailed  description  of  the  LINEAR  observing  strategy  is  given 
in  Appendix  A. 

2.3.  Photometry  and  Astrometry  Using  SExtractor 

The  LINEAR  project  has  been  archiving  all  of  its  image 
data  since  2002  December  through  the  present.  The  data  used 
for  this  work  encompass  all  imagery  collected  from  2002 
December  through  2008  March.  After  selecting  the  LINEAR 
imagery  data  that  overlap  the  main  SDSS  survey  region  (i.e., 
galactic  latitudes  greater  than  30°)  and  SDSS  stripe  82  region 
(20h32m  <  ay20oo.o  <  04h00m,  -126  <  8jmm)  <  +E26, 
~280  deg2),  the  data  comprise  about  1.8  million  images.  The 
bias  and  flat-field  corrections  were  applied  to  these  images 
following  standard  methodology. 

The  extraction  of  sources  from  LINEAR  images  and  the  cre¬ 
ation  of  source  catalogs  were  done  using  the  SExtractor  software 
(Bertin  &  Arnouts  1996;  see  the  SExtractor  user’s  manual6  for 
more  details  on  the  source  detection  and  measurement  process). 
The  preliminary  astrometry  for  extracted  sources  was  obtained 
by  converting  centroid  pixel  coordinates  into  equatorial  J2000.0 
right  ascension  (R.A.)  and  declination  (decl.)  by  application  of 
the  gnomonic  projection  and  a  priori  telescope  pointing  infor¬ 
mation.  This  projection  does  not  correct  for  telescope  pointing 
errors  on  the  order  of  several  pixels,  telescope  optical  distortion 
on  the  order  of  a  few  pixels  near  the  corners  of  the  image,  or 
differential  refraction  on  the  order  of  a  few  pixels  at  typical  air- 
masses.  When  compared  to  the  USNO-B  catalog  (Monet  et  al. 
2003)  astrometry,  the  average  uncertainty  in  the  preliminary 
astrometry  is  about  in  R.A.  and  decl.  coordinates.  There¬ 
fore,  the  astrometry  reported  by  SExtractor  must  be  recalibrated 
before  matching  extracted  LINEAR  sources  against  other  cat¬ 
alogs.  The  recalibration  of  preliminary  astrometry  is  described 
in  Section  3. 

The  preliminary  fixed-aperture  photometry  for  each  source 
is  computed  by  summing  the  pixels  within  a  5  pixel  (Ml") 
radius  circular  aperture  and  subtracting  the  background  estimate 
for  each  pixel  computed  in  the  background  estimate  step.  The 
uncertainty  on  the  computed  flux  is  computed  assuming  Poisson 
statistics  in  detected  electrons  with  a  gain  conversion  factor  of 
2.5e~  ADU-1.  The  pixels  with  counts  above  17,000  were  not 
used  as  the  CCDs  begin  to  experience  nonlinear  response  above 


6  Available  at  https://www.astromatic.net/pubsvn/software/sextractor/ 
trunk/doc/ sextractor.pdf. 


this  level  (corresponding  to  42,000  electrons),  and  they  saturate 
at  full- well  capacity  of  75,000  electrons. 

3.  RECALIBRATION  OF  LINEAR 
ASTROMETRY  AND  PHOTOMETRY 

In  this  section,  we  describe  the  astrometric  and  photometric 
recalibration  of  LINEAR  data  and  quantify  the  quality  of 
recalibrated  data.  The  recalibration  and  data  quality  assessment 
are  done  on  a  field-by-field  basis,  where  a  field  is  defined 
as  part  of  a  LINEAR  image  observed  by  a  single  CCD  cell 
(there  are  four  cells  in  the  small  setup  and  eight  cells  in 
the  LI  and  L2  setups).  The  only  exception  to  this  is  the 
astrometric  recalibration,  which  is  done  on  a  image-by-image 
basis.  However,  the  quality  of  astrometric  recalibration  is  still 
assessed  on  a  field-by-field  basis. 

3.1.  Astrometric  Recalibration 

We  recalibrate  the  preliminary  astrometry  using  the  Astrome- 
try.net  software  (Lang  et  al.  2010).  Astrometry.net  uses  the  x  and 
y  (hereafter,  pixel)  positions  of  sources  extracted  by  SExtractor 
and  looks  for  “skymarks”  (sets  of  four  stars)  that  correspond  to 
pre-computed  “skymarks”  found  in  the  “clean”  USNO-B  cat¬ 
alog  (Barron  et  al.  2008).  Once  a  sky  mark  is  found  the  code 
obtains  an  initial  astrometric  solution  (a  transformation  from 
the  pixel  coordinates  into  the  J2000.0  Equatorial  system)  in  the 
form  of  a  third-degree  polynomial  in  v  and  y  coordinates.  Using 
this  initial  astrometric  solution,  the  code  looks  for  additional 
reference  USNO-B  stars  in  the  image,  and  further  refines  the 
initial  astrometric  solution  when  it  finds  them. 

The  quality  of  recalibrated  LINEAR  astrometry  is  estimated 
by  comparing  the  LINEAR  positions  with  reference  USNO-B 
positions  (Figure  3).  We  find  that  the  recalibrated  LINEAR 
single-epoch  coordinates  have  an  rms  scatter  of  MX' 6  and  are 
offset  <(X3  with  respect  to  the  USNO-B  coordinates.  These 
systematic  offsets  are  independent  of  pixel  positions  or  the 
position  of  the  image  on  the  sky  and  are  subtracted  from 
recalibrated  LINEAR  coordinates.  Using  a  subset  of  data,  we 
have  found  that  the  offsets  decrease  if  Astrometry.net  is  allowed 
to  search  for  even  more  reference  USNO-B  stars  when  refining 
the  initial  astrometric  solution. 

To  find  fields  with  bad  astrometry,  we  calculate  the  median 
and  rms  scatter  of  astrometric  residuals  (difference  between 
LINEAR  and  USNO-B  positions)  in  R.A.  and  decl.  for  each 
field,  and  tag  fields  as  bad  if  the  median  or  rms  scatter  in 
either  coordinate  (R.A.  or  decl.)  is  greater  than  2".  We  find  that 
about  7%  of  all  fields  have  bad  astrometry,  with  the  majority  of 
bad  fields  (~90%)  originating  from  the  L2  setup,  and  remove 
them  from  further  processing.  Visual  inspection  of  images 
containing  bad  fields  reveal  that  bad  fields  lack  several  pixel 
columns,  which  were  most  likely  dropped  due  to  cell  readout 
problems. 

In  order  to  recalibrate  the  LINEAR  photometry,  we  posi¬ 
tionally  matched  the  LINEAR  data  set  to  the  SDSS  Data  Re¬ 
lease  7  (SDSS  DR7;  Abazajian  et  al.  2009)  imaging  catalog. 
The  SDSS  DR7  imaging  catalog  provides  homogeneous  and 
deep  (r  <  22.5)  l%-2%  accurate  photometry  in  five  band- 
passes  ( u ,  g ,  r,  i,  and  z)  of  more  than  11,000  deg2  of  the  North 
Galactic  Cap,  and  three  stripes  in  the  South  Galactic  Cap  total¬ 
ing  740  deg2.  The  SDSS  positions  are  accurate  to  better  than 
Qf'.l  per  coordinate  (rms)  for  sources  with  r  <  20.5  (Pier  et  al. 
2003),  and  the  morphological  information  from  the  images  al¬ 
lows  reliable  star-galaxy  separation  to  r  ~  21.5  (Lupton  et  al. 
2002). 
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R-A-iinear  -  R-A.usnob  [arcsec] 


DecIinear  -  Decusnob  [arcsec] 

Figure  3.  Comparison  of  USNO-B  and  recalibrated  LINEAR  R.A.  (top)  and 
decl.  single-epoch  coordinates  (bottom)  for  the  small  (solid  line),  LI  (dotted 
line),  and  L2  (dashed  line)  setups.  The  width  of  the  distributions  (determined 
from  the  interquartile  range)  puts  the  astrometric  precision  in  LINEAR  single¬ 
epoch  coordinates  at  ~(X'6.  The  systematic  offsets  in  LINEAR  coordinates  are 
<(X'3  and  have  been  subsequently  corrected  for.  The  dotted  line  at  0"  was  added 
to  guide  the  eye. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

We  match  LINEAR  sources  to  SDSS  DR7  non-saturated, 
“primary”  objects7  using  a  2"  radius  and  find  ~4.8  billion 
matched  pairs  (there  are  24  million  unique  SDSS  sources 
included  in  these  pairs).  There  are  ~1.6  billion  LINEAR 
sources  without  an  SDSS  match  (“orphans”),  and  they  include 
observations  associated  with  stars  saturated  in  SDSS  (r  <  13, 
~60%  of  orphans),  image  artifacts  (e.g.,  due  to  charge  bleeding 
and  diffraction  spikes),  moving  objects  such  as  asteroids,  and 
other  transients. 

3.2.  Photometric  Recalibration 
We  model  the  recalibrated  LINEAR  magnitudes  as 

^linear  =  O^instr  +  C(g  ~  /)  +  Z(x,  y),  (1) 

where  a  measures  nonlinearity  of  SExtractor- supplied  instru¬ 
mental  magnitudes  minstr  =  —  2.5  log(counts),  C  is  a  color- 
dependent  term  that  accounts  for  the  per-exposure  change  in 
the  effective  central  wavelength  of  the  LINEAR  bandpass,  and 
Z(x,y)  is  our  model  for  the  field- specific  magnitude  zero-point 
dependent  on  the  x-y  position  in  the  image,  defined  as 

Z(x,  y)  =  Zo  +  Z\x  +  Z^y  +  F(x,  y),  (2) 

where  Zo,  Z\,  and  Z2  vary  from  field  to  field,  and  F(x,y)  is 
the  “super  flat  field”  (fixed  for  each  observational  setup).  Note 


7  See  http://cas.sdss. org/dr6/en/help/docs/algorithm.asp?key=flags. 


Figure  4.  Color  term  between  the  zero-point-corrected  LINEAR  magnitude, 
Winstr  +  Zo,  and  SDSS  r-band  magnitude  for  the  small  (dots),  LI  (solid  squares), 
and  L2  (open  triangles)  setups.  The  solid  line  shows  the  best  fit  for  all  three 
setups.  The  rms  scatter  between  the  best  fit  for  all  three  setups  and  fits  to 
individual  setups  is  smaller  than  0.01  mag.  The  SDSS  g  —  i  color  is  not  corrected 
for  ISM  extinction. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

that  the  color-independent  airmass  term  (gray  extinction)  is 
implicitly  included  in  Z0. 

To  recalibrate  the  LINEAR  photometry,  we  need  a  catalog 
of  photometric  standard  stars  in  the  LINEAR  photometric 
(unfiltered)  system.  Since  such  a  catalog  does  not  exist,  we  use 
the  SDSS  imaging  catalog  and  utilize  SDSS  gri  photometry  to 
derive  magnitudes  in  the  LINEAR  photometric  system. 

3.2.1.  LINEAR  Magnitudes  Synthesized  from  SDSS  Photometry 

We  define  LINEAR  magnitudes  synthesized  from  the  SDSS 
photometry  as 

Wsdss  =  r  +  f{g-  i),  (3) 

where  f(g  —  i)  is  some  function  of  g  —  i  and  SDSS  gri 
magnitudes  are  not  corrected  for  the  ISM  extinction. 

We  find  f(g  —  i )  iteratively  by  first  adopting  a  =  1.0  in 
Equation  (1),  F(x,  y)  =  0  in  Equation  (2),  and  f(g  —  i)  =  0 
in  Equation  (3).  For  each  field,  we  determine  Zo,  Zi,  and  Z2  by 
minimizing 

X  =  ^[(^linear  ^sdss)/°rlinear]  5  (4) 

where  a\mciir  is  the  magnitude  error  supplied  by  the  SExtractor 
(errors  in  msdss  are  assumed  to  be  much  smaller  than  errors  in 

^linear)- 

The  LINEAR  observations  of  calibration  stars  must  be 
matched  within  2"  to  isolated,8  unresolved,  SDSS  objects 
brighter  than  r  —  17,  with  0  <  g  —  i  <  2.6.  Having  Z(x,y) 
for  each  field,  we  calculate  mm str  +  Z(x,  y)  —  r  residuals  for  all 
calibration  stars,  bin  the  values  in  g  —  i  bins,  and  plot  the  median 
residuals  in  Figure  4  for  the  small,  LI,  and  L2  setups.  We  plot 
medians  for  each  setup  separately  because  each  setup  represents 
a  separate  optical,  and  therefore  photometric,  system. 

Even  though  the  three  setups  are  different,  their  residuals 
depend  similarly  on  g  —  i  (rms  scatter  around  the  mean  is 
0.01  mag),  as  indicated  by  medians  plotted  in  Figure  4.  We 
therefore  fit  a  fourth-degree  polynomial  to  all  medians  and 
obtain 

fig  -  i)  =  0.0574  +  0.004(g  -  i)  -  0.056(g  -  if 

+  0.052(g  -if  -  0.0262 (g  -if.  (5) 


SDSS  CHILD  flag  set  to  0. 
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Figure  5.  Distribution  of  raiinear  —  ^sdss  residuals  binned  in  the  x-y  plane  for  the  small  (top),  LI  (middle),  and  L2  (bottom)  setups.  The  panels  on  the  left  show  median 
values  in  each  bin,  and  the  panels  on  the  right  show  rms  scatter  in  a  bin.  The  values  are  color-coded  according  to  legends,  with  values  outside  the  range  saturating.  The 
patterns  seen  on  the  left  are  magnitude  zero-point  variations  that  were  not  taken  out  when  LINEAR  images  were  flat  fielded.  Therefore,  the  maps  on  the  left  represent 
“super  flat  fields”  that  need  to  be  subtracted  during  recalibration. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


The  transformation  from  SDSS  gri  to  LINEAR  magnitudes,  as 
defined  by  Equation  (3),  is  then 

msdss  =  r  +  0.0574  +  0.004(g  —  i)  —  0.056(g  —  if 

+  0.052(g  -  if  -  0.0262(g  -  if.  (6) 

This  equation  effectively  turns  the  SDSS  imaging  catalog  into  a 
catalog  of  LINEAR  photometric  standard  stars. 

3.2.2.  Super  Flat  Fields  and  Nonlinearity  of  Instrumental  Magnitudes 

We  find  the  super  flat  field,  E(v,y),  iteratively  by  first  adopting 
a  =  1.0  in  Equation  (1)  and  F(x,y)  =  0  in  Equation  (2). 
For  each  field,  we  determine  Zo,  Zi,  Z2,  and  C  by  minimizing 
Equation  (4)  using  calibration  stars  brighter  than  msdss  =  17. 
The  super  flat-field  correction,  F(x,y ),  is  then  obtained  by  finding 


the  median  of  raiinear  —  ^sdss  residuals  binned  in  x-y  pixels, 
and  is  shown  as  a  color-coded  map  in  Figure  5  (left)  for 
each  observational  setup.  Figure  6  (left)  shows  the  residuals 
after  the  subtraction  of  Figure  5  (left)  maps.  On  average,  the 
median  residuals  are  now  smaller  than  0.005  mag,  compared  to 
~0. 02-0. 03  mag  before  the  correction  was  made. 

To  verify  our  assumption  of  a  =  1.0  (the  linearity  of 
instrumental  magnitudes),  we  adopt  a  =  1.0  in  Equation  (1), 
and  determine  Z0,  Z\,  Z2,  and  C  by  minimizing  Equation  (4) 
using  calibration  stars  brighter  than  msdss  =  17.  We  bin 
^linear  —  ^sdss  residuals  in  msdss  bins,  and  plot  the  medians 
as  a  function  of  msdss  in  Figure  7.  A  strong  dependence  of 
residuals  on  msdss  would  indicate  nonlinear  ( a  /  1.0)  behavior 
of  instrumental  magnitudes.  With  the  exception  of  a  small 
nonlinearity  (<0.02  mag)  in  the  LI  setup  for  msdss  <  15, 
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Figure  6.  Similar  to  Figure  5,  but  obtained  after  subtraction  of  “super  flat  fields.”  The  median  values  are  color  coded  on  a  smaller  range  than  in  Figure  5  (0.02  mag 
vs.  0.05  mag).  On  average,  the  medians  are  now  smaller  than  0.005  mag,  and  the  rms  scatter  is  smaller  than  0.05  mag.  The  higher  rms  scatter  indicates  cells  with 
sometimes  problematic  readouts. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


the  residuals  do  not  depend  on  magnitude,  indicating  that  the 
a  =  1.0  assumption  is  valid. 

To  estimate  the  magnitude  of  correction  introduced  by  in¬ 
cluding  the  color-dependent  term  C  (accounts  for  the  change  in 
LINEAR  bandpass  due  to  varying  airmass)  in  Equation  (1),  we 
measure  the  rms  scatter  of  C(g  —  i)  values  obtained  from  all  cali¬ 
bration  stars  and  fields.  We  find  the  rms  scatter  to  be  ~0.02  mag, 
indicating  that  by  correcting  for  this  term  the  systematic  uncer¬ 
tainty  in  LINEAR  photometry  is  reduced  by  ~0.02  mag.  On 
the  other  hand,  the  inclusion  of  the  Zo  term  is  much  more  im¬ 
portant  as  it  reduces  the  systematic  uncertainty  in  LINEAR 
photometry  by  ~0.3  mag  (the  rms  scatter  of  Z0  values  is 
~0.3  mag). 

Once  the  color-dependent  term,  C,  and  the  zero  point,  Z(v,y), 
are  obtained  for  a  field,  the  LINEAR  photometry  for  that  field  is 
recalibrated  using  Equation  (1).  For  “orphan”  sources  (LINEAR 


sources  not  associated  with  SDSS  objects)  the  color-dependent 
term  is  set  to  C  =  0,  as  g  —  i  colors  are  not  available  for  these 
sources.  Therefore,  the  systematic  uncertainty  in  recalibrated 
photometry  is  slightly  higher  for  orphans  (by  ~0.02  mag). 

3.2.3.  Checks  on  Recalibrated  Photometry 

To  quantify  the  quality  of  recalibration  on  a  field-to-field 
basis,  we  calculate  miinear  —  wsdss  residuals  for  calibration  stars, 
and  find  their  median  and  rms  scatter  (determined  from  the 
interquartile  range)  for  each  field.  The  median  of  (miinear  —  msdss) 
residuals  per  field  is  then  an  estimate  of  the  zero-point  error. 
The  distribution  of  these  median  offsets,  shown  in  Figure  8,  is 
about  0.01  mag  wide  and  demonstrates  the  overall  quality  of 
recalibrated  photometry.  About  10%  of  fields  are  of  low  quality 
(rms  scatter  in  mnmar  —  msdss  >0.1  mag),  usually  due  to  highly 
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^sdss 

Figure  7.  Medians  of  m\inear  —  ms dss  residuals  in  msdSs  bins  for  the  small  (dots), 
LI  (crosses),  and  L2  setups  (open  triangles).  Apart  from  slight  non-linearity 
(<0.02  mag)  in  the  LI  setup  at  bright  magnitudes,  recalibrated  LINEAR 
magnitudes  are  linear  with  magnitudes  synthesized  from  SDSS  gri  photometry. 
(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


Figure  8.  Distribution  of  zero-point  offsets  between  SDSS  and  LINEAR  systems 
for  the  small  (solid  line),  LI  (dotted  line),  and  L2  (dashed  line)  setups,  where 
the  offset  is  measured  for  each  field  as  the  median  of  miinear  —  ^sdss  residuals. 
The  distribution  is  about  0.01  mag  wide,  implying  that  on  average  the  LINEAR 
recalibration  is  about  0.01  mag  offset  from  the  SDSS  photometric  system. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


variable  cloud  coverage.  Sources  detected  in  such  fields  are 
removed  from  further  consideration. 

Another  useful  statistic  is  the  rms  scatter  of  /  values, 
(^linear  —  ^sdss) /^linear,  in  ^sdss  bins.  This  rms  scatter  should 
be  equal  to  1  and  should  not  depend  on  msdss  if  the  SExtractor- 
supplied  magnitude  errors  (^linear)  are  correct  Gaussian  errors. 
As  shown  in  Figure  9,  the  rms  scatter  is  greater  than  1  indicating 
underestimated  aimear  errors  by  up  to  a  factor  of  ~2. 

We  correct  the  underestimated  magnitude  errors  by  multi¬ 
plying  them  by  1.4,  1.85,  and  1.65  for  the  small,  LI,  and  L2 
setups,  respectively.  In  the  case  of  the  small  setup,  the  factor  of 
1 .4  is  only  the  average  value  as  SExtractor-supplied  magnitude 
errors  for  that  setup  seem  to  depend  on  magnitude.  Therefore, 
the  magnitude  errors  for  the  small  setup  may  be  slightly  over¬ 
estimated  or  underestimated  depending  on  the  brightness  of  the 
source.  This  is  not  a  major  issue  as  only  a  small  fraction  of  all 
LINEAR  sources  (~0.1%)  originate  from  the  small  setup. 

The  reason  for  this  underestimation  of  SExtractor-supplied 
magnitude  errors  is  not  clear.  Non-inclusion  of  systematic  errors 
(e.g.,  due  to  flat  fielding)  in  aiinear  cannot  explain  this  behavior 
for  LI  and  L2  setups.  If  additive  systematic  errors  were  the 
reason  for  this  underestimation,  then  (miinear  —  msdss) Minear 
should  decrease  with  magnitude  because  aiinear  is  expected  to 
increase  with  magnitude  (see  Figure  13).  Yet,  no  such  decrease 
is  seen  in  Figure  9  for  LI  and  L2  setups. 


Figure  9.  rms  scatter  of  (m\\near  —  msdss) /^linear  values  in  msdss  bins  for  the 
small  (dots),  LI  (crosses),  and  L2  (open  triangles)  setups.  The  rms  scatter 
values  greater  than  1  indicate  underestimated  LINEAR  magnitude  errors.  The 
pipeline-supplied  errors  for  the  small,  LI,  and  L2  setups  are  underestimated  by 
about  40%,  85%,  and  65%,  respectively. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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Figure  10.  Median  number  of  observations  per  object  as  a  function  of  equatorial 
J2000.0  right  ascension  and  declination  coordinates.  The  values  are  color-coded 
according  to  the  legend,  with  values  outside  the  range  saturating.  The  dashed 
lines  show  ±10°  of  the  ecliptic  plane  (solid  line)  and  the  Galactic  plane  is  shown 
as  a  dotted  line.  The  average  number  of  observations  per  object  within  ±10°  of 
the  ecliptic  plane  is  ~460,  and  ~200  elsewhere. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

4.  CHARACTERIZATION  OF  THE  RECALIBRATED 
LINEAR  DATA  SET 

There  are  about  5  billion  sources  in  the  recalibrated  LINEAR 
data  set.  These  sources  were  grouped  into  clusters  using  an 
implementation  of  the  OPTICS  algorithm  (Ankerst  et  al.  1999) 
and  a  clustering  radius  of  2".  There  are  about  25  million  clusters 
with  more  than  15  observations.  SDSS  provides  morphological 
classification  for  24  million  clusters  (7  million  are  resolved 
and  17  million  are  unresolved  by  SDSS),  and  for  1  million 
of  these  clusters  there  is  no  SDSS  classification  (they  are 
slightly  outside  the  SDSS  footprint).  Clusters  with  15  or  more 
sources  are  hereafter  labeled  as  LINEAR  objects  or  time  series. 
The  median  number  of  epochs  (observations)  per  object  as 
a  function  position  and  magnitude  is  shown  in  Figures  10 
and  11.  The  average  number  of  epochs  for  objects  brighter 
than  18  mag  within  ±10°  of  the  ecliptic  plane  is  ~460,  and 
~200  elsewhere. 

Clusters  with  less  than  15  sources  are  not  tagged  as  objects. 
We  have  statistically  analyzed  properties  of  385  million  sources 
(~8%  of  all  LINEAR  sources)  assigned  to  such  clusters  and 
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Figure  11.  Median  number  of  “good”  observations  (observations  with  flags  = 
0)  as  a  function  of  msdSs-  The  median  number  of  “good”  observations  is  ~200 
for  msdss  <18,  after  which  it  decreases  precipitously  (the  5 o  detection  limit  is 
at  msdss  ~  18).  About  10%  of  detections  have  flags  >0. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


Xpdf 


Figure  12.  Cumulative  distribution  of  x^f  values  for  a  bright  (14.5  < 

linear)  <  17,  solid  line)  and  faint  (17  <  (miinear)  <  18,  dashed  line)  sample  of 
LINEAR  objects.  The  median  of  distributions  is  ~1,  indicating  that,  on  average, 
the  uncertainties  in  LINEAR  magnitudes  are  well  determined.  The  high-end  tail 
of  distributions  (Xpdf  >  3)  is  populated  by  truly  and  spuriously  variable  objects 
(most  likely  dominated  by  the  latter). 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

have  found  that  most  of  them  are  very  faint  (miinear  >  19.5), 
have  very  uncertain  magnitudes  (crimes  >0.4  mag),  and  have 
high  ellipticity  (suggestive  of  cosmic  rays  and  image  artifacts). 
Such  clusters  and  their  observations  are  not  included  in  our 
database.  This  also  means  that  the  LINEAR  database  does  not 
include  observations  of  fast  moving  objects  (i.e.,  objects  moving 
faster  than  ~1"  day-1). 

In  order  to  enable  the  discovery  of  variable  objects  and  to 
further  characterize  the  LINEAR  data  set,  we  have  computed 
various  low-order  statistics  for  time  series,  such  as  the  median 
recalibrated  magnitude  ((miinear))*  rms  scatter  (a),  x2  Per 
degree  of  freedom  (Xpdf),  skewness  (yi),  and  kurtosis  ( yi )•  The 

cumulative  distributions  of  Xpdf  for  LINEAR  objects  with  more 
than  30  observations  are  shown  in  Figure  12. 

The  median  of  the  Xpdf  distribution  for  bright  LINEAR  objects 
(14.5  <  (miinear)  <  17)  is  about  1,  indicating  that,  on  average, 
the  uncertainties  in  LINEAR  magnitudes  are  well  determined. 
For  fainter  LINEAR  objects  (17  <  (miinear)  <  18),  the  median 
of  the  Xpdf  distribution  is  slightly  lower  (~0.7)  indicating 
slightly  overestimated  photometric  errors. 

The  high-end  tail  (x2df  >  3)  of  the  cumulative  Xpdf  dis¬ 
tribution  suggests  that  ~5%  of  bright  LINEAR  objects  ex¬ 
hibit  deviations  from  the  mean  that  are  inconsistent  with  the 
Gaussian  distribution  of  errors.  These  objects  could  be  variable 


Figure  13.  Median  photometric  error,  X,  as  a  function  of  magnitude.  The  median 
photometric  error  is  estimated  as  the  median  of  o  values  in  a  magnitude  bin 
(dots),  where  o  is  the  rms  scatter  of  a  time  series.  The  analytic  description  of 
£("flinear)  is  given  by  Equation  (7)  (solid  line),  and  was  obtained  by  fitting  a 
fourth-degree  polynomial  to  medians.  The  median  photometric  error  increases 
rapidly  from  about  0.03  mag  at  15  mag,  where  systematic  errors  dominate,  to 
0.2  mag  at  18  mag,  where  the  errors  are  dominated  by  photon  statistics. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

or  could  have  unreliable  (non-Gaussian)  errors.  While  it  is  not 
impossible  that  some  fraction  of  LINEAR  objects  are  variable 
(~1%  of  sources  brighter  than  r  =  20  are  variable;  see  Sesar 
et  al.  2006,  and  references  therein),  a  fraction  of  ~5%  seems 
unlikely.  A  more  likely  explanation  is  spurious  variability  in¬ 
duced  by  non-Gaussian  errors,  caused  by  blending  of  sources 
and  unaccounted  instrumental  defects  (e.g.,  bad  pixels). 

To  estimate  the  dependence  of  the  median  photometric  error 
(X)  on  magnitude,  we  follow  Sesar  et  al.  (2007)  and  calculate 
median  a  values  in  (miinear)  bins.  The  dependence  of  median 
o  on  (miinear),  as  shown  in  Figure  13,  can  be  modeled  as  a 
fourth-degree  polynomial: 

n.x)  =  122.995  -  30.346.x  +2.81 83x2  -  0.116893.x3 

+  0.0018295.x4,  (7) 

where  x  =  miinear.  The  smallest  photometric  errors  are  at 
^linear  ~  15  (~0.03  mag)  and  are  limited  by  systematic  errors 
in  calibration.  The  increase  in  photometric  error  at  brighter 
magnitudes  is  probably  due  to  saturation  issues  in  LINEAR. 
The  photometric  error  starts  to  increase  rapidly  toward  fainter 
magnitudes  due  to  photon  noise  and  reaches  ~0.2  mag  at 
^linear  ~  18  (the  adopted  faint  limit). 

In  addition  to  statistics  calculated  on  light  curves,  we  also 
calculated  time-averaged  positions  (median  and  mean  R.A.  and 
deck).  To  estimate  the  reliability  of  time-averaged  LINEAR  po¬ 
sitions,  we  compared  them  to  SDSS  positions  in  ~1  deg2  pixels 
on  the  sky  and  calculated  the  median  difference  in  positions 
for  each  pixel.  We  find  that,  on  average,  the  time-averaged 
LINEAR  positions  are  within  073  of  SDSS  positions,  with  no 
systematic  trends  across  the  sky.  The  largest  median  differ¬ 
ence  (for  a  ~1  deg2  pixel)  is  ~1"5.  While  this  precision  is 
sufficient  for  most  applications,  we  still  recommend  that  SDSS 
positions  are  used  whenever  possible  to  avoid  patches  where 
time-averaged  LINEAR  positions  may  be  less  reliable. 

Since  SDSS  is  about  4  mag  deeper  than  LINEAR,  it  can 
be  used  to  estimate  the  completeness  of  the  LINEAR  ob¬ 
ject  catalog  (i.e.,  LINEAR  object  catalog  is  a  subset  of  the 
SDSS  object  catalog).  To  estimate  the  completeness  of  the 
LINEAR  object  catalog,  in  Figure  14  we  plot  the  fraction  of 
SDSS  objects  matched  to  LINEAR  objects  as  a  function  of  msdss 
(LINEAR  magnitude  synthesized  from  SDSS  photometry;  see 
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Figure  14.  Completeness  of  the  LINEAR  catalog  as  a  function  of  msdss 
magnitude  (LINEAR  magnitude  synthesized  from  SDSS  photometry,  see 
Equation  (6)).  The  completeness  is  estimated  as  a  fraction  of  SDSS  objects 
matched  to  LINEAR  objects  within  2" .  We  find  the  object  catalog  to  be  more 
than  90%  complete  for  msdss  <19. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 


Equation  (6)).  For  this  purpose,  we  use  SDSS  and  LINEAR 
objects  with  0  <  g  —  i  <2.5  (color  range  where  Equation  (6) 
is  valid),  170°  <  R.A.  <  220°,  and  15°  <  decl.  <  65°.  As 
shown  in  Figure  14,  the  LINEAR  object  catalog  is  more  than 
90%  complete  for  msdss  <19.  This  estimate  may  be  a  bit  mis¬ 
leading  as  the  criterion  for  object  creation  is  having  at  least  15 
detections.  Even  faint  objects  (18  <  ms dss  <  19.5)  are  likely  to 
be  detected  at  least  15  times  out  of  an  average  of  200  observa¬ 
tions  (see  Figure  1 1),  if  the  observing  conditions  are  better  than 
average  (e.g,  due  to  better  seeing  and  darker  skies). 


5.  PRELIMINARY  ANALYSIS  OF  VARIABLE  OBJECTS 
SELECTED  FROM  THE  LINEAR  DATA  SET 

We  have  examined  a  small  subset  of  objects  with  Xpdf  >  3 
(~5000)  and  identified  a  subset  of  periodically  variable  stars. 
We  confirmed  expectation  from  the  preceding  section  that  only 
a  few  percent  of  the  Xpdf  >  3  subsample  are  convincing 
cases  of  periodic  variability.  As  an  illustration  of  the  quality 
of  LINEAR  light  curves  at  the  bright  ((miinear)  ~  15)  and 
faint  ((miinear)  ~  17)  ends,  we  show  their  period-folded  light 
curves  in  Figure  15.  The  main  strength  of  the  LINEAR  survey 
is  easily  discernible  in  these  plots:  even  at  the  faint  end,  the 
light  curve  features  are  well  defined,  thanks  to  dense  LINEAR 
sampling. 

To  exploit  the  full  potential  of  the  LINEAR  photometric 
database,  an  automated  selection  and  classification  of  variable 
objects  is  desirable.  This  desideratum  is  not  easy  to  accomplish 
because  the  fraction  of  spurious  objects  with  Xpdf  >  3  is  an  order 
of  magnitude  larger  than  the  fraction  of  truly  variable  objects. 
In  order  to  enable  the  deployment  of  automated  selection 
and  classification  methods,  we  have  undertaken  an  extensive 
program  of  visual  light  curve  classification.  Details  of  this 
program  are  described  in  a  companion  paper  (L.  Palaversa 
et  al.  201 1,  in  preparation,  hereafter  Paper  II)  and  here  we  only 
provide  a  brief  description  of  the  variable  candidate  selection 
and  their  distribution  in  various  color-color  diagrams. 

In  the  first  step,  we  define  a  flux-limited  sample  of  objects 
with  14.5  <  (miinear)  <  17  and  select  about  200,000  light  curves 
with  Xpdf  >  3  and  a  >  0. 1  mag  (a  is  the  rms  scatter  in  the  light 
curve).  We  limit  classification  to  objects  exhibiting  periodic 
variability  and  use  phased  light  curves  for  visual  inspection. 


Phase 


Phase 


Figure  15.  Examples  of  period-folded  light  curves  from  the  LINEAR  data  set.  The  left  panels  show  bright  ((munear)  ~  15)  and  right  panels  show  faint  objects 
((^linear)  ~  17).  The  two  top  light  curves  are  RR  Lyrae  stars,  and  other  are  eclipsing  binary  systems.  The  number  in  the  lower  left  corner  in  each  panel  indicates  the 
rms  scatter  around  the  cubic  spline  fit  to  the  binned  mean  values  of  the  light  curve  using  20  points  (dashed  line).  The  light  curve  features  are  well  defined,  thanks  to 
dense  LINEAR  sampling,  even  for  faint  stars. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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<  -0.7  -0.55  -0.4  -0.25  >  -0.1 

log(period/day) 


Figure  16.  Distribution  of  log(period)  for  ~7000  LINEAR  variables  in  the  SDSS 
u  —  g  vs.  g  —  r  color-color  diagram.  The  sample  is  limited  to  14.5  <  r  <  17 
and  variability  was  confirmed  by  visual  classification  of  light  curves.  Sources 
are  binned  in  0.05  mag  wide  color  bins,  and  the  median  values  are  color- 
coded  according  to  the  legend  on  the  right  (values  outside  the  range  saturate 
in  blue  or  red).  Contours  outline  the  overall  distribution  of  SDSS  stars  (with 
14.5  <  r  <  17)  at  1%,  10%,  and  50%  of  the  peak  density. 

(A  color  version  of  this  figure  is  available  in  the  online  journal.) 

The  periods  are  determined  using  the  Supersmoother  algorithm 
(Reimann  1994).  In  total,  about  7000  LINEAR  objects  (3.5% 
out  of  200,000)  were  visually  confirmed  as  periodic  variables, 
with  about  two- thirds  classified  as  RR  Lyrae  type  ab  stars  and 
eclipsing  binary  stars  of  all  types.  The  remaining  1/3  include 
SX  Phoenicis  and  Delta  Scuti  stars,  long-period  variables,  as 
well  as  possible  RR  Lyrae  type  ab  and  eclipsing  binary  stars 
with  less  certain  classification.  The  distributions  of  classified 
objects  in  the  SDSS  u  —  g  versus  g  —  r  color-color  diagram  is 
shown  in  Figure  16. 

The  differences  between  the  distributions  of  all  sources  and 
those  of  the  confirmed  variable  subsample  demonstrate  that 
the  latter  are  robustly  classified  (i.e.,  the  classification  is  not 
random).  The  most  obvious  difference  is  a  much  higher  fraction 
of  RR  Lyrae  stars  in  the  classified  sample  (u  —  g  ~  1.15, 
g  —  r  <  0.3).  The  eclipsing  binary  systems  of  all  types  dominate 
the  visually  classified  sample  in  the  main  stellar  locus  region. 
There  seems  to  be  a  hint  for  the  median  period  of  eclipsing 
binaries  becoming  longer  for  bluer  stars.  The  distribution  of 
eclipsing  stars  is  offset  from  the  distribution  of  all  stars  in  the 
direction  perpendicular  to  the  stellar  locus  (offset  to  upper  left). 
There  is  also  a  lack  of  variable  (eclipsing)  systems  of  the  late-K 
and  early-M  spectral  type.  These  results  are  discussed  in  more 
detail  in  Paper  II. 

6.  SUMMARY 

In  this  work,  we  have  described  astrometric  and  photometric 
recalibration  of  data  obtained  by  the  asteroid  survey  LINEAR. 
The  public  access  to  the  recalibrated  LINEAR  data  set  will  be 
provided  through  the  SkyDOT  Web  site  (http://skydot.lanl.gov). 
A  description  of  various  cataloged  parameters  that  are  available 
in  this  database  is  given  in  Appendix  B. 

The  preliminary  astrometry  produced  by  the  survey  has  been 
recalibrated  using  the  Astrometry.net  software  and  USNO-B 


catalog  as  the  reference  astrometric  catalog.  The  precision 
(rms  scatter)  of  single-epoch-recalibrated  LINEAR  coordinates 
is  ~0T6  with  respect  to  the  reference  (USNO-B)  astrometric 
catalog.  The  precision  of  time-averaged  LINEAR  coordinates 
is  ~( y/3  with  respect  to  SDSS  positions. 

In  order  to  recalibrate  the  LINEAR  photometry,  we  position- 
ally  matched  the  LINEAR  data  set  to  the  SDSS  DR7  imag¬ 
ing  catalog  and  used  SDSS  as  a  catalog  of  photometric  stan¬ 
dard  stars.  LINEAR  instrumental  magnitudes  were  corrected 
for  changes  in  the  LINEAR  bandpass  (due  to  varying  airmass) 
and  for  zero-point  variations  across  the  CCD  (due  to  imperfect 
flat  fielding).  These  corrected  LINEAR  magnitudes  were  then 
offset  on  a  per-LINEAR  field  basis  to  match  the  SDSS  mag¬ 
nitude  zero-point  system.  In  general,  the  recalibrated  LINEAR 
magnitudes  (raiinear)  behave  linearly,  with  the  exception  of  a 
small  nonlinearity  (<0.02  mag)  in  the  LI  setup  at  magnitudes 
brighter  than  15  mag. 

We  find  the  average  uncertainty  in  recalibrated  LINEAR 
magnitudes  to  be  ~0.03  mag  at  the  bright  end  (miinear  ~  15,  the 
systematic  uncertainty)  and  ~0.2  mag  at  the  faint  end  (miinear  ~ 
18,  the  adopted  faint  limit  for  the  survey).  The  SExtractor- 
supplied  magnitude  errors  were  found  to  be  underestimated  by 
40%-85%,  independent  of  magnitude,  and  were  subsequently 
corrected  for  this  behavior  using  multiplicative  factors. 

There  are  about  5  billion  photometric  measurements  in  the 
recalibrated  LINEAR  data  set.  The  sources  detected  in  each 
observation  were  grouped  into  clusters  (time  series  or  objects) 
using  a  clustering  radius  of  2".  There  are  about  25  million 
clusters  with  at  least  15  observations  (7  million  resolved  and 
17  million  unresolved  by  SDSS,  the  rest  without  morphological 
classification).  The  median  number  of  epochs  (observations) 
per  object  brighter  than  18  mag  and  within  ±10°  of  the  ecliptic 
plane  is  ~460,  and  ~200  elsewhere.  The  completeness  of  the 
LINEAR  object  catalog  (relative  to  the  SDSS  imaging  catalog) 
was  estimated  to  be  better  than  90%  for  sources  brighter  than 
19  mag. 

LINEAR  data  set  represents  a  major  new  resource  for 
studying  the  variability  of  faint  optical  sources.  With  an  average 
of  200  observations  per  object,  LINEAR  data  provide  time 
domain  information  for  the  brightest  four  magnitudes  of  the 
SDSS  survey.  At  the  same  time,  LINEAR  significantly  extends 
the  deepest  similar  wide-area  variability  survey,  the  Northern 
Sky  Variability  Survey,  by  3  mag. 

In  order  to  help  automated  discovery  of  variable  objects, 
we  have  computed  various  low-order  statistics  for  time  series 
(median  recalibrated  magnitude,  rms  scatter,  x2  per  degree 
of  freedom,  skewness,  and  kurtosis).  The  Xpdf  distribution  for 
bright  LINEAR  objects  (brighter  than  17  mag)  is  centered  at 
~1  (indicating  that,  on  average,  the  uncertainties  in  LINEAR 
magnitudes  are  well  determined),  and  has  a  long  high-end 
(X^jf  >  3)  tail.  About  5%  of  LINEAR  objects  have  x^f  >  3;  in 
a  sample  with  Gaussian  errors,  less  than  0.01%  of  objects  are 
expected  to  have  x^f  >  3.  This  high-end  tail  is  populated  by 
truly  variable  objects,  but  is  most  likely  dominated  by  spuriously 
variable  objects  with  unreliable  (non-Gaussian)  errors  (caused 
by,  for  example,  unaccounted  instrumental  defects  such  as  bad 
pixels  or  by  blended  sources). 

Since  the  fraction  of  spurious  variable  objects  is  at  least  a 
few  times  greater  than  the  fraction  of  truly  variable  objects, 
we  have  undertaken  an  extensive  program  of  visual  light  curve 
classification  in  order  to  enable  future  deployment  of  automated 
methods.  Details  of  this  program  are  described  in  a  companion 
paper  (L.  Palaversa  et  al.  201 1,  in  preparation). 
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APPENDIX  A 

LINEAR  OBSERVING  STRATEGY 

To  search  for  asteroids,  the  LINEAR  program  attempts  to 
cover  the  entire  portion  of  sky  available  from  the  site  each  lunar 
month  (lunation),  with  repeat  searches  within  10°  to  15°  of 
the  ecliptic  plane.  The  program  schedules  searches  on  both 
telescopes  an  average  of  25  nights  per  lunation,  taking  off 
4-5  nights  around  the  full  moon.  After  accounting  for  weather 
and  equipment  problems,  about  20  observing  nights  are  achieved 
per  lunation  with  each  telescope.  To  minimize  stray  light,  the 
survey  fields  are  kept  about  20°  away  from  the  moon  during  the 
week  before  and  after  the  full  moon. 

To  efficiently  cover  the  available  sky  at  minimal  airmass,  each 
telescope  is  assigned  a  strip  of  sky  that  is  aligned  with  either  the 
ecliptic  or  the  equator  and  is  typically  3-4  fields  of  view  (4°-5°) 
wide  in  declination  or  ecliptic  latitude,  and  150  fields  of  view 
(240°)  long  in  right  ascension  or  ecliptic  longitude.  The  search 
box  is  divided  into  5-7  regions,  each  containing  about  100  fields. 
The  regions  are  searched  starting  with  the  westernmost  region 
at  the  beginning  of  astronomical  twilight  in  the  evening,  and 
ending  with  the  easternmost  region  at  the  end  of  astronomical 
twilight  in  the  morning.  All  of  the  fields  in  a  single  region  are 
imaged  in  sequence,  and  then  the  sequence  is  repeated  in  that 
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region  four  more  times.  The  five  images  of  a  given  field  are 
combined  into  a  “frameset”  which  forms  the  input  to  a  moving 
target  detection  algorithm  that  distinguishes  asteroids  from  stars 
and  noise  based  on  their  motion.  This  process  of  cycling  through 
the  fields  in  a  region  five  times  produces  a  temporal  separation 
of  15-20  minutes  between  successive  images  of  the  same  patch 
of  sky. 

On  subsequent  nights,  the  entire  search  area  for  each  telescope 
is  generally  stepped  north  or  south  by  the  height  of  the  region 
in  order  to  cover  the  entire  sky,  though  gaps  are  occasionally 
unavoidable  due  to  missed  nights.  Over  the  course  of  one 
lunation  the  entire  available  sky  is  generally  searched  once, 
and  areas  within  about  10°  of  the  ecliptic  are  searched  two  or 
three  times.  Because  the  monthly  search  regions  are  quite  wide 
in  right  ascension,  there  is  significant  overlap  in  coverage  area 
on  adjacent  lunar  months.  Therefore,  a  given  patch  of  sky  can 
be  observed  up  to  eight  months  in  a  row  before  it  passes  into  the 
daytime  portion  of  the  sky  for  four  months.  After  a  full  year  the 
search  covers  the  entire  sky  north  of  declination  —30°  except 
for  a  small  circle  north  of  declination  82°  where  the  telescope 
has  poor  settling  ability.  LINEAR’s  pattern  of  sky  coverage 
produces  repeat  photometry  data  for  stars  on  several  timescales 
ranging  from  the  15-20  minute  interval  between  images  within  a 
frameset,  to  a  few  days  between  repeat  visits  during  one  lunation, 
to  the  month-long  timescale  between  lunar  months,  to  the 
yearly. 

APPENDIX  B 

LINEAR  DATABASE  SCHEMA 

The  LINEAR  database  hosted  at  SkyDOT  contains  two  tables, 
“Object”  and  “Source.”  The  “Source”  table  contains  photome¬ 
try  and  astrometry  of  individual  LINEAR  detections  (sources), 
while  the  “Object”  table  references  time  series  obtained  by  clus¬ 
tering  sources  within  2".  In  addition,  the  “Object”  table  contains 
various  low-level  statistics  calculated  on  time  series  and  is  cross- 
matched  to  external  catalogs  such  as  SDSS  and  Two  Micron  All 
Sky  Survey  (Skrutskie  et  al.  2006).  A  detailed  description  of 
various  parameters  provided  by  these  tables  is  given  in  Tables  1 
and  2,  and  a  screenshot  showing  the  SkyDOT  interface  is  shown 
in  Figure  17. 


Table  1 

LINEAR  “Source”  Table  Schema 


Field 

Description 

sourcelD 

A  unique  number  identifying  this  source 

ra 

Equatorial  J2000.0  right  ascension 

decl 

Equatorial  J2000.0  declination 

MJD 

Modified  Julian  Date 

mag 

Recalibrated  LINEAR  magnitude 

magErr 

Uncertainty  in  the  recalibrated  LINEAR  magnitude 

calib 

1  if  the  source  was  used  for  calibration,  0  otherwise 

ellip 

“Ellipticity”  parameter  measured  by  SExtractor 

ximg 

Position  of  the  source  along  CCD’s  x-axis 

yimg 

Position  of  the  source  along  CCD’s  y-axis 

FWHM 

FWHM  measured  by  SExtractor 

flagsa 

SExtractor  internal  flags 

objectID 

A  number  associating  this  source  with  an  object  in  the  “Object”  table 

Note.  a  See  p.  28  in  the  SExtractor  manual  (https://www.astromatic.net/pubsvn/software/sextractor/ 
trunk/doc/ sextractor.pdf) . 
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Table  2 

LINEAR  “Object”  Table  Schema 


Field 


Description 


objectID 

ra 

decl 

ra_std 

decl_std 

ra_median 

decl_median 

nptsTotal 

nptsGood 

mag_mean 

mag_wmean 

mag_median 

mag_std 

mag_wstd 

mag_skew 

mag_kurtosis 

mag_chi2dof 

mag_rchi2dofa 

theta_sdssb 

sdssID 

sdss_ra 

sdss_dec 

objtype 

rExt 

uMod 

gMod 

rMod 

iMod 

zMod 

uErr 

gErr 

rErr 

iErr 

zErr 

ug 

gr 

ri 

iz 

gi 

pmLc 

pmBc 

pmErrc 

isolated 

sdss_satur 

sdss_flags 

theta_2massd 

J 

J  En¬ 
id 

HErr 

K 

KErr 

jh 

hk 

jk 


A  unique  number  identifying  this  object 
Mean  equatorial  J2000.0  right  ascension 
Mean  equatorial  J2000.0  declination 
rms  scatter  in  right  ascension  (in  units  of  arcsec) 
rms  scatter  in  declination  (in  units  of  arcsec) 

Median  equatorial  J2000.0  right  ascension 
Median  equatorial  J2000.0  declination 
Number  of  detections 
Number  of  detections  with  flag==0 
Mean  recalibrated  LINEAR  magnitude 

Weighted  recalibrated  LINEAR  magnitude  (by  magnitude  errors) 

Median  recalibrated  LINEAR  magnitude 
rms  scatter  of  recalibrated  LINEAR  magnitudes 

Weighted  rms  scatter  of  recalibrated  LINEAR  magnitudes  (by  magnitude  errors) 

Skewness  of  recalibrated  LINEAR  magnitudes 
Kurtosis  of  recalibrated  LINEAR  magnitudes 
X2  per  degree  of  freedom  around  the  mean  magnitude 
Robust  x2  Per  degree  of  freedom  around  the  mean  magnitude 

Distance  in  arcsec  from  the  SDSS  source 
SDSS  Data  Release  7  objID  parameter 
SDSS  J2000.0  right  ascension 
SDSS  J2000.0  declination 

6  if  unresolved  in  SDSS  imaging  (i.e.,  a  point  source),  3  if  resolved 

Extinction  in  the  SDSS  r-band 

SDSS  w-band  model  magnitude 

SDSS  g-band  model  magnitude 

SDSS  r-band  model  magnitude 

SDSS  /-band  model  magnitude 

SDSS  z-band  model  magnitude 

SDSS  w-band  model  magnitude  error 

SDSS  g-band  model  magnitude  error 

SDSS  r-band  model  magnitude  error 

SDSS  /-band  model  magnitude  error 

SDSS  z-band  model  magnitude  error 

SDSS  u  —  g  color  (corrected  for  ISM  extinction,  u  —  g  =  (wMod  —  1.87rExt)  —  (gMod  —  1.38rExt)) 
SDSS  g  —  r  color  (corrected  for  ISM  extinction,  g  —  r  =  (gMod  —  1.38rExt)  —  (rMod  —  l.OOrExt)) 
SDSS  r  —  /  color  (corrected  for  ISM  extinction,  r  —  i  =  (rMod  —  l.OOrExt)  —  (/Mod  —  0.76rExt)) 
SDSS  /  —  z  color  (corrected  for  ISM  extinction,  iz  =  (/Mod  —  0.76rExt)  —  (zMod  —  0.54rExt)) 
SDSS  g  —  i  color  (corrected  for  ISM  extinction,  g  —  /  =  (gMod  —  1.38rExt)  —  (/Mod  —  0.76rExt)) 
Proper  motion  in  the  Galactic  latitude  direction  (in  units  of  mas  yr_l) 

Proper  motion  in  the  Galactic  longitude  direction  (in  units  of  mas  yr-1) 

Proper  motion  error  (in  units  of  mas  yr-1) 

1  if  SDSS  CHILD  flag  is  set  to  0  (isolated  ==  0  may  indicate  a  close  neighbor) 

1  if  SDSS  SATURATED  flag  is  set  to  1 
SDSS  flags  encoded  as  a  64-bit  integer  number 

Distance  in  arcsec  from  the  2MASS  source 
2MASS  7-band  magnitude 
2MASS  7-band  magnitude  error 
2MASS  H- band  magnitude 
2MASS  H- band  magnitude  error 
2MASS  K- band  magnitude 
2MASS  A'-band  magnitude  error 

2MASS  7  — //color  (corrected  for  extinction,  j  —  h  =  (7  —  0.327rExt)  —  (H  —  0.209rExt)) 

2MASS  H  —  K  color  (corrected  for  extinction,  h  —  k  =  (H  —  0.209rExt)  —  (K  —  0.133rExt)) 
2MASS  7—  K color  (corrected  for  extinction,  (7  —  0.327rExt)  —  (K  —  0.133rExt)) 


Notes. 

a  Calculated  after  eliminating  10%  of  highest  and  10%  of  lowest  magnitudes. 

b  This  and  other  SDSS  parameters  set  to  NULL  if  no  match  within  2"  and  not  in  SDSS  footprint,  set  to  100  if  in  SDSS  footprint 
but  no  match  within  2" . 

c  Obtained  from  the  Munn  et  al.  (2004)  SDSS-USNOB  proper  motion  catalog. 

d  This  and  other  2MASS  parameters  set  to  NULL  if  no  match  within  2".  Following  Covey  et  al.  (2007),  we  only  match  to  2MASS 
sources  with  the  following  2MASS  flags:  rd^ag  =  222,  bl^ag  =  111,  and  ccflag  =  0. 
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Now  Viewing:  user_sql_query 


SQL  query 


select  •  from  source  where  "obiectlD"  =  1030; 

vvwvvvwv\ 


Output  Online  Table  C 


!  Reset 
Run  Query 

Your  query:  "select  ♦  from  source  where  "objectID"  =  1030;  " 
Matching  records 


item 

sourcelD 

ra 

decl 

MJD 

mag 

magErr 

calib 

ellip 

1 

1634337 

115.843476 

65.736482 

52636.3064595 

14.07 

0.01 

0 

0.025 

2 

3209021 

115.843379 

65.736453 

52636.3498227 

14.105 

0.01 

0 

0.005 

3 

6385157 

115.843423 

65.736468 

52636.3209077 

14.077 

0.01 

0 

0.035 

4 

6508400 

115.843364 

65.736455 

52636.3353559 

14.074 

0.01 

0 

0.08 

5 

6815434 

115.843422 

65.736473 

52636.2917809 

14.069 

0.01 

0 

0.058 

6 

10495189 

115.843327 

65.736382 

52621.3346657 

14.066 

0.008 

0 

0.037 

1  1 

in 

J 

io 

'7 

FI 

Page  1 

|  Of  19  f  >i| 

Displaying  1  to  10  of  186  items 
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Figure  17.  Screenshot  showing  the  search  interface  for  the  LINEAR  database. 
(A  color  version  of  this  figure  is  available  in  the  online  journal.) 
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